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We derive an explicit expression for the Fourier-Laplace transform of the two-point distribution 
function p{x\,t\; X2,t2) of a continuous time random walk (CTRW), thus generalizing the result of 
Montroll and Weiss for the single point distribution function p{xi,ti). The multi-point distribution 
function has a structure of a convolution of the Montroll- Weiss CTRW and the aging CTRW single 
point distribution functions. The correlation function {x{ti)x{t2)) for the biased CTRW process is 
found. The random walk foundation of the multi-time-space fractional diffusion equation [Baule 
and Friedrich [Europhysics Letters 77 10002 (2007)] is investigated using the unbiased CTRW in 
the continuum limit. 

PACS numbers: 05.40.Fb, 02.50.-r 



I. INTRODUCTION 

The continuous time random walk (CTRW) introduced 
in [l[ more than forty years ago is a successful model for 
normal and anomalous diffusion in a variety of physical 
systems 0, H, 0, H, 0] ■ A new splash of interest in this old 
model was caused by the fact that it provides a firm sta- 
tistical foundation of the fractional Fokker-Planck equa- 
tion 0, [Sj , and is a simple model for the investigation of 
such intriguing phenomena of non-equilibrium statisti- 
cal physics as weak ergodicity breaking and statistical 
aging [13, [m, [H, [11. The importance of CTRW as a 
minimal model of non-Markovian behavior is connected 
with its semi-Markovian (renewal) nature, which allows 
in many cases for an exact probabilistic description of the 
process of interest. Unlike Markovian diffusion processes, 
which are fully characterized by their transition probabil- 
ities, non-Markovian CTRW requires for the full hierar- 
chy of multi-point distribution functions for its complete 
characterization [H, [H, [H, [l3, [il| . Multi-point distri- 
bution functions and correlation functions are necessary 
tools to distinguish between CTRW stochastic dynam- 
ics from other non Markovian processes, e.g. fractional 
Brownian motion, as appearing in linear fracton mod- 
els or in generahzed Rouse models, see e.g. [l9|. The 
simplest experimental way of discriminating these types 
of anomalous processes is based on different behavior of 
their correlation functions, as investigated e.g. in recent 
experiments on anomalous diffusion in single molecules 
[20! . [2ll . [2^ . Hence obviously it is a worthy goal to 
consider higher-order joint probability density functions 
(PDFs) and correlation functions of the CTRW. 



In the present article we concentrate on decoupled 
CTRWs in the sub-diffusive and spatially homogeneous 
limit and obtain an exact expression for the two time 
characteristic function of the CTRW process in Laplace 
space. This main result is a generalization of the 
MontroU-Weiss equation [ij which gives the character- 
istic function of the single time PDF. We use our result 
to investigate the validity of a multi-point fractional ki- 
netic equation of Baule and Friedrich ^Ig] . We moreover 
obtain the two-point correlation function {xi{ti)x2{t2)) 
in a biased CTRW. We note that even this simple correla- 
tion function cannot be found from the Green's function 
p{xi,ti) of the CTRW, since the process is non Marko- 
vian. 



II. DEFINITIONS AND NOTATION 

We consider a standard CTRW model in one dimen- 
sion with a walker starting at the origin at time t — Q 
[E [1, [3, 13, [1, [1 ■ Waiting times between jumps are inde- 
pendent identically distributed (IID) random variables 
with a common PDF il}{t). After waiting the particle 
makes a jump and the process is then renewed. The 
jump lengths are also IID random variables, with a PDF 
f{Sx). The waiting times and the lengths of jumps are 
mutually independent. 

Let p{xi,ti;x2,t2) = p{xi,ti;x2,t2\0,0) be the joint 
PDF of finding a walker at xi at time ti > and at X2 at 
time t2 > 0. We concentrate first on the corresponding 
multi-point characteristic function, i.e. on the double 
Fourier double Laplace transform of p{xi, ti; X2, 12): 



p(fci, Ai; fc2, A2) 



dxi 



dx-2 



dti 



dt2 e 



p{xi,ti;X2,t2). 
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All over the article, the original functions and their trans- 
forms can be distinguished on the ground of their vari- 
ables (xi, ti; X2, and (fti, Ai; fc2, A2) for originals and 
transforms, respectively. The same holds for other func- 
tions encountered in the text. 

Let us now separate in Eq.([T]) the contributions corre- 
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sponding to the cases ti < t2 and ti > t2 and rewrite the 
integral as a sum of the two terms: 

p(fci,Ai;fc2,A2) = p<(fci,Ai;fc2,A2) + p>(A:i, Ai; /c2, A2) 

(2) 

with 



p<(fci, Ai;A;2, A2) = 
p>(fci, Ai;A:2, A2) = 



dxi 
dxi 



00 

dh 
Jt 

poo I'll 

dx2e'^^''^+'^^''^ / dh I dt2e 
-00 Jo Jo 



p{xi,ti;X2,t2) 

^^*'~^'''p{xi,h;x2,t2). 



(3) 
(4) 



For ti < t2 we define the elapsed time t — t2 — ti > and the corresponding displacement A = 12 — Xi. Similarly, 
for ti > t2 we take t = ti — t2 and A = xi — X2. For t2 > ti we define the joint PDF g{xi,ti; A, t) that the particle 
is at xi at time ti and then experiences a displacement of size A during the time interval r. Similar joint PDF for 
t2 < ti is denoted with g(x2,t2; A, r). 

To obtain p{ki, Ai; ^2, A2) we use g{xi,ti; A, r) and g{x2, ^2; A, r) in the corresponding terms of Eq. ([3]) and Eq. 
([4]). Changing in Eq. ([3]) the variables according to ^2 = + t and X2 — xi + A one gets 



/OO pQC pOO />OC 

dxi / dA / dh / dre-^i*i-^^(*i+")+^'=^"i+*^-^("^+^)g(xi,ii;A,T) 
-00 J-oc Jo Jo 



g{ki + fc2,Ai + A2;fc2,A2). 



(5) 



A similar change of variables is made in Eq. ^ giving us the corresponding displacements 



p(A:i, Ai;A;2, A2) = g(fci + /c2, Ai + A2; fci, Ai) 

-I- 5(A:i +/c2,Ai + A2;fc2,A2). (6) 

Thus, our next task is to find the characteristic func- 
tion g. We consider ti < t2 and soon concentrate on 
g{ki, Xi; k, s) which is the double Laplace and double 
Fourier transform of g{xi,ti] A,t) according to the rule 
xi ki, ti ^ Xi, A —>■ k and t —> s. Let ni be the ran- 
dom number of jumps made by a walker during the time 
interval (0, ti), n2 the number of jumps made in the time 
interval (ti, ti+r), and Pni,n2 (^ii ''') the joint probability 
of these random variables. Since waiting times and jump 
lengths are independent, we may write 

00 

g{xi,ti;A,T)= ^ Pn,,n2{ti,T)p{xi;ni)p{A\xi;n2) 

ni ,712—0 

(7) 

where p{xi; ni) is the probability density to find a parti- 
cle at xi after ni steps and p{A\xi; 712) is the probability 
density of the displacement A of a particle starting at 
xi after ri2 steps. Eq. ([7]) is the key to all our further 
considerations, and it shows that the problem can be di- 
vided into three: calculation of i^ni.n2 (^ii ''') which is not 
trivial, and the calculation of p{xi;ni) and p{A\xi;n2). 
The latter two describe simple random walks in discrete 
time, moreover, since the jump lengths are independent 
p{A\xi;n2) does not depend directly on xi. In this case 



ni~\-n2 



xi=^6xi, A= ^ 6xi, 



(8) 



z— ni + l 



are sums of IID random variables, and their characteristic 
functions are: 



(g^«i^i) = I p(xi;ni)e'^''='dxi = (e*'=i*^)"i = /"^(fci) 

(9) 



and 



p(A|xi;n2)e*'-^dA = r^(fc) (10) 



where /(fc) ~ {exp{ik5x)) is the one step characteristic 
function i.e. the Fourier transform of f{5x). Passing to 
the Fourier ~ Laplace representation we hence get 

00 

g(fci,Ai;fc,s)= Pn,..2{>^i,s)n{k^)r{k). 

,712—0 

(11) 

where P„j^„2(Ai,s) is the double Laplace transform of 
^Tii, 712(^1 1 ''') in its temporal variables. We thus see that 
this Laplace-transform PTii,n2(-^ii ^) of the probability of 
the number of steps during the corresponding time inter- 
vals plays the key role in our theory and we turn now to 
determining it. 
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III. STATISTICS OF NUMBERS OF STEPS 

-P"l,"2(-^ll ^) 

Let us now consider the set of jump times as a point 
process on the time axis and let ti^t2, ■ ■ -U ■ ■ ■ denote 
the corresponding points: ti is the instant of time when 
the first jump was made, t2 is the time when the sec- 
ond jump was made etc. As mentioned, according to the 
CTRW model the waiting times ti, t2 ~ ti, — <2 etc 
are IID random variables with the common PDF tpit)- 
The realizations of the process with ni jumps up to time 
ti and jumps between ti and ti + t are those that 
satisfy t„j < ti < tn^+i and t„i+„2 < ti + t < t„^+n2+i 
respectively. We introduce the indicator function I(x) 
of a logical variable x so that I{x) = 1 if the condition 
X in the brackets holds {x is true) and I{x) — other- 
wise. As usual the probability for x to be true is then 
the mean {I{x)) of I{x) over realizations. For the case 
under consideration 

Pn„nAtl,T)= (12) 

ni -{-712 

The argument of the second indicator function can be 
rewritten as {tm+n2 - h < t < t„j+„2+i - h)- The 
double Laplace transform of P„j.„2 (ii, i") is 



-Pni,ri2('^lj s) 



/ dre~^" / dtie"^^*^P„„„,(ti,r) 
JO Jo 



dhe'^^'H (tn, <ti< tn, + l) X 



J dre "'^ I {tn^+n2 - tl < T < tni+n2 + l - tl)"^ {^3) 

where we used the linearity of the Laplace transform to 
interchange the sequence of integration and ensemble av- 
eraging. Note that only the second indicator function 
contains r as a variable and the second integral in Eq. ([Te 



S.I. = / dre '^'^ I [tni+n2 - h < T < tni+n2 + l - h) 

Jo 

is rather trivial and it exhibits three behaviors: (i) S.L = 

for t„i+„2+i - <i < 0, (m) 



S.L 



1 -e 



-s(t„j + „2 + l-tl) 



for i„i+„2 < ^1 < ini+na+i and {Hi) 

"^1 + 112 + 1-*! 



S.L 



s(t™i+„2-*l) _ g-s(t..l + ..2 + l-*l) 



for tl < tni+7i2- We note that condition (ii) and the 
condition tm < ti < tm+i can hold simultaneously only 
when 712=0, while condition (iii) and the condition tm < 
tl < tni+i can only hold simultaneously if n2 7^ 0. Using 
these behaviors we now get from Eq. ([Tc 



-Pni,o(Ai, s) 



tni + l 



dtie-''^U(tm <ti< tm 



_ -s{tm + i-ti) 



1 - e 



+1 



(14) 



for 772 = and 



/ rt,^l+,^2 _ _ g-s(t'n+"2^*l) _ g ^ « ( i + 712 + 1 - * 1 ) 

fni,n2(Ai,s) = (J dtie~'^^H {tm < tl < tm+i) 

/ p-st„i+„2 _ g-St„j_|_„2_|_i p-(Ai-s)t„j _ g-(Ai-s)t„j + i \ 



(15) 



for 7i2 7^ 0. 



Since waiting times are IID random variables, one has 
(e-*"i^i) =7^"MAi) (16) 

and 

-st»ii + ig-(Ai-s)t„j \ _ / -Ait„jg-s(t„j + i-t„j) 



(17) 



V'"^(Ai)V(s). 



Similar expressions hold also for other terms in Eqs. (|T^ 
and Here -0 (Ai) and V'(s) are Laplace transforms of 
the waiting time PDF '0(^1 ) and V'(t) respectively. Using 
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these expressions we get: 

S 

and 

-Pni,n2 ('^Ij = 



l-^(Ai) ^(s)-^(Ai) 



Ai 



(18) 
(19) 



(Ai - s) 



for n2 > 1. Note that Eqs. (|18|19p give the proper nor- 
mahzation since EJ^^^q I^i^=o -Pai,s("i, "2) = l/ (Ais). 

We now consider limiting behaviors of Eqs. (|18ll9p . 
The probabihty of making no steps in the time interval 
{ti, ti + t) is given in double Laplace representation by 



00 



E 

ni=0 



X P„i,o(Al,s) = — 
— ' sAi 



V'(s)-V(Ai 



s(Ai-s) [l-V^(Ai)]^ 



(20) 



which was obtained previously [23]. Let t/ be the time 
between ti and the first jump event after ti: tf = 
tni+i — ti. The random variable tf is sometimes called 
the forward recurrence time. Let its PDF be given by 
h(tf]ti) depending on ti as a parameter. In double 
Laplace representation ti — > Ai and tf s one finds 



h{s\\i) 



jjjs) - ^(Ai) 
(Ai-s) [l-^(Ai)]^ 



(21) 



as follows from Eq. (pOI by noting that the probability 
of making no jump in the time interval (ti,ii + r) is 
1 — h{tf]ti)c\tf . The probability of making n2 > 1 
jumps in (ti,ti + r) is 



P„i,„2 (Ai, s) = h(s\ Ai) 

rii=0 



l-^(s) 



^"^-i(s). (22) 



In the space of originals this equation corresponds to a 
convolution of the PDF of forward recurrence time tf 
with the PDFs of the the following n2 — 1 waiting times; 
the factor [1 — ip{s)]/s comes from the probability of not 
jumping between the last event in the sequence and the 
end of observation at <2- 



IV. TWO-POINT CHARACTERISTIC 
FUNCTIONS 

We are now able to find the characteristic function 
g(fci,Ai;fc,s), Eq. pT|. using Eqs. (|18|19p : 



g(fci, Ai;fc,s) 



l-^(Ai) V(s)-^(Ai) 



1 



Ais s (Ai — s) 

f{k) [l-V(s)] [^(s)-^(Ai)] 



l-^(Ai)/(fci) 



1 



1 



s(Ai 
1 fe(s;Ai) 
Ais s 



1 



l"^(Ai)/(fci) l-V(s)/(fc) 
V'(Ai) , l-^(Ai) /i(s;Ai)/(fc) l-^(s) 



l-V(Ai)/(fci) l-^(Ai)/(Ai) 



l-^(s)/(fc)- 



(23) 



This equation can be written in a more transparent 
way. First we recall the MontroU-Weiss equation. For 
a CTRW starting at time t = 0, the single point PDF 
Pyiwi^jt) of the particle being at site x at time t, is 
given in Laplace t —>■ X Fourier x —>■ k space in terms of 
the MontroU- Weiss equation [l| 



i^MW [^(A),/(A)] = 



1 - ^(A) 



1 



A l-/(fc)V(A)' 



(24) 



The MontroU- Weiss equation explicitly assumes that the 
waiting time for the first step has the same PDF as all 
further waiting times. On the other hand, one may con- 
sider situations where the waiting time PDF for the first 
step, '01 (t) differs from the PDFs of all other waiting 
times 'ip{t) [2^. The single point PDF describing this 
more general process [Tol, nil called aging random walk 
[2^ is denoted by PARwix,t). In the Laplace- Fourier 



representation a; — > fc, t ^ A one finds 

PARw[^i(A),^A(A),/(fc)] = (25) 
l-V'i(A) , Vi(A)/(fc) l-^(A) 



A 



l-/(A:)V'(Ai) A 



The aging random walk reduces to the MontroU- Weiss 
CTRW if 01 (t) = 0(t). Using Eqs. ^ and Eq.^ one 
can rewrite Eq. ((23)) as 



g{ki,Xi;k,s)= (26) 
Pmw [V- (Ai) , / (fci)] Parw [Ai/i(s; Ai), V (s) , / (fc)] . 

We see that the solution for g corresponds to a convolu- 
tion of two PDFs, the one of the MontroU- Weiss CTRW 
and the one of the aging CTRW with the first waiting 
time PDF formally put to iPi{t) = dh{T,t)/dt. We note 
that the fact that the final characteristic function is a 
convolution and not a simple product of Pmw (2^1, ii) and 
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Parw(^i T]ti), as found for simple Markovian diffusion, 
has to do with the correlations between xi and X2 which 
in turn is related to the correlation between ni and n2- 
These arise through subtle correlations between the num- 



ber of steps ni and the forward recurrence time. 

After getting g{ki, Ai; fc, s) we can turn to the original 
characteristic function. Using Eqs. ^ [^5]) we find 



Ai; fc2, A2) = + A2), /(fci + ^2)] Parw [(Ai + A2) /i(Ai + A2; A,), (A,) , / (h)] . (27) 

I 



One can check that if fci = or ^2 = (i.e. integrating 
the overall distribution over xi or X2 respectively) we 
recover the MontroU- Weiss equation (p4)) for a one-point 
characteristic function, for example 



p{ki = 0, Ai; A;2, A2) = 



1 



1-V'(A2) 



A1A2 l-^(A2)/(fc2) 



(28) 



as it should. 

We now investigate the continuum limit of our main 
equation (|7f|) . corresponding to long ti and ^2, using the 
standard long wave length (small k) and small frequency 
approximation Q. We consider first the non-biased ran- 
dom walks with a finite second moment of jump lengths 



J 



((5x^), which means that for small k 

{Sx^)P 



For Ai 



/(fc)^l 



^(Ai) ~ 1- A(Ai)^ 



(29) 



(30) 



where < a < 1 and yl > 0. For the special case a — 
1, yl is the mean waiting time. This case corresponds 
to asymptotically normal diffusion. If a < 1 the mean 
time between jumps diverges, which leads to anomalous 
behaviors. In this limit 



p(fci, Ai; A;2, A2) 



(A1 + A2) 



Q-l 



(Xi + X2T + D^\ki + k2\ 



E 



(Ai + A2)(A,)"-A, (Ai + A2)^ 



A1A2 (Ai-hA2) 



+ (Ai + A2) 



(Ai+A2)"-(A0" (Ai)°-i 
Aj(Ai + A2)" (A,)" + 2?„(fc,)2 



(31) 



where 1 = 2 and 2 = 1 and Da = {Sx'^) / {2A) is the fractional diffusion constant ,8]. Previously Baule and Friedrich 
[1^ wrote a multi-point fractional diffusion equation for p{xi,ti; X2,t2), whose solution in Laplace- Fourier space is 
exactly Eq.(l31]). ^ 



V. CORRELATION FUNCTION FOR BIASED 
CTRW 



We now consider the simplest correlation function 



(xi(Ai)a;2(A2)) 



d d 
dki dk2 



p(fci, Ai; ^2, A2)|fci=fc2=o. 



(32) 

For a biased CTRW with finite variance of jump lengths, 
the small k expansion reads 



/(fc) ~ l + i{5x)k 



(33) 



where {5x) is the mean step length. Using Eqs. 
(I27I32I33I) 



(a;i(Ai)a;2(A2)) = 



A1A2 [l-^(Ai+A2) 

2{5x)^ 
A1A2 [l-^(Ai+A2) 



(34) 



{5xf 



1 - 



fe(Ai + A2;A,) 



E 



Now we pass to the small Ai and A2 limit and consider 
the scaling limit of large ii and t2 when their ratio is 
arbitrary, using Eqs. pOI34p 



(a;i(Ai)a;2(A2)) 



AAiA2(Ai +A2)" 
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jSx)' l/(Ai)" + l/iX,Y 
AiA2(Ai + A2)" 



(35) 



Laplace transform of Eq. (j35p is found using the ap- 
proach discussed in Appendix B of Ref. 17]. For t2 > ti 
one has: 



{Xi{ti)x2{t2)) 



A r(l + a; 



(36) 



r(i + 2a) r(i 



F a 



t2 



Taking ti > t2 corresponds to simple interchange of the 
arguments. Here F(a,b]c; z) is a hypergeometric func- 
tion [13] ■ Let us check limiting behaviors of the correla- 
tion function Eq. ([55)1 . For an unbiased process {6x) ~ 
the first term on the right hand side of equation ([36]) is 
the only non- vanishing term 



{Xi{ti)x2(t2)) 



A r(l + a) 



(37) 



which represents the dispersion of the walker's positions 
at time ti. On the other hand, for the biased CTRW 
{Sx) 7^ the leading behavior is given by the second line 
in Eq. (|36p and we may neglect the first term. Several 
limiting situations can be used as benchmarks. For t2 ^ 
h one has F{a, -a; a -I- 1; x) ~ 1 — a'^x/ (1 -|- a) -|- 0{x'^) 
so that 



with 



{Xi{h)x2{t2)) {Xl{h)){x2{t2)) 



A r(l-l-a) 
A r{l + a)' 



(38) 



(39) 



which proves the decoupling of correlations for t2 3> ti. 
In the opposite limit of t2 ti one uses Eq. (15.3.6) of 
Ref.([27|) a.nd F{a,b;c;0) = 1 to get i^(a, -a; a + 1; 1) = 
[r(a + l)]2/r(2a + 1) and to obtain 



lim {xi{ti)x2{t2)) 



{6x^) {hr , 2{5x)\t,f 



A r(l + a) r(l + 2a)^2' 

(40) 

which is the mean square displacement in the biased 
CTRW d^, as expected. 

Another limit is the Markovian case a = 1, for which 
F(l,— l;2;a:) = 1 — x/2. In this case the decoupling, 



Eq. ([38[) . is valid at all times. The fact that for a < 1 
Eq. ([55|) holds only for t2 ^ ii, i.e. the existence of non- 
trivial correlations between xi and X2 , has to do with the 
correlations between the number of steps before and after 
the first observation time ti which we discussed in Sec. 
IV. Thus, if relatively few jumps take place during the 
time interval (0,ti), i.e. ni <C {ni), the typical displace- 
ment xi is inevitably small, and then the particle is likely 
to be effectively trapped at its position at ti for a very 
long time which is of the order of ti. In this case also 
the forward recurrence time is long in statistical sense. 
This implies that also n2 is going to be relatively small, 
since the particle will likely wait for a long time for its 
first step after ti, which leads to small absolute values of 
X2 as well. Hence correlations for a < 1 are built even 
when ti and t2 are very long. 



VI. DISCUSSION 



The MontroU- Weiss equation ([M]) expresses the char- 
acteristic function of the CTRW in terms of Laplace and 
Fourier transforms of the PDFs of the waiting times and 
jump lengths. Similarly, Eq. ([27)1 gives the two di- 
mensional characteristic function of the CTRW process. 
From this equation we may derive two dimensional cor- 
relation functions for the CTRW process, for example we 
considered the biased CTRW. We showed that the two 
dimensional characteristic function depends on the prob- 
ability of ni renewals in (0, ti) and 712 renewals in (ii, 12) 
(for t2 > ti) and that these numbers of steps are corre- 
lated. For characteristic functions of order N higher than 
two one would have to calculate renewal statistics in 
intervals. In principle this calculation can be performed 
using the same technique we used here for example to 
calculate Pni,„2,n3(ii,i2,t3) etc. 

The two dimensional characteristic function is shown 
to be related to the MontroU- Weiss and aging CTRW 
single point characteristic functions. Thus even though 
the process is non-Markovian information on one dimen- 
sional characteristic functions is sufficient to find the two 
dimensional characteristic function. This simplification 
is obviously related to the renewal property of the un- 
derlying random walk. Finally, starting with the CTRW 
model we derived the solution of the multi-point frac- 
tional diffusion equation [Tsj . in Fourier-Laplace space, 
thus giving further justification for this new equation. 
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